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Abstract. We have searched for solar axions or similar particles that couple to two 
photons by using the CERN Axion Solar Telescope (CAST) setup with improved 
conditions in all detectors. From the absence of excess X-rays when the magnet 
was pointing to the Sun, we set an upper limit on the axion-photon coupling of 
< 8.8 X 10-"GeV"^ at 95% CL for ma < 0.02 eV. This result is the best 
experimental limit over a broad range of axion masses and for ma < 0.02 eV also 
supersedes the previous limit derived from energy-loss arguments on globular-cluster 
stars. 
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1. Introduction 

The Peccei-Quinn (PQ) mechanism [1-3] remains perhaps the most compelling 
explanation of the CP problem of strong interactions. The existence of a new chiral 
f/(l) symmetry that is spontaneously broken at some large energy scale /a would allow 
for the dynamical restoration of the CP symmetry in QCD. An inevitable consequence of 
this mechanism is the existence of axions, the Nambu-Goldstone bosons of U (1)pq [4, 5]. 
Axions are pseudoscalar particles with properties closely related to those of neutral 
pions. In particular, their mass and interaction strengths can be obtained by a simple 
rescaling from the corresponding 7r° properties (see equations (3) and (4) below). 
Unsuccessful experimental searches [6] and astrophysical limits [7] imply that axions, if 
they exist, must be very light and very weakly interacting, yet they are not necessarily 
harmless because they are a leading candidate for the cold dark matter of the universe [8- 
11]. Moreover, several approaches to a unified description of dark energy and dark 
matter involving pseudo Nambu-Goldstone bosons have been put forward recently [12- 
14]. 

Experimental strategies to find "invisible axions" or more general "axion-like 
particles" (ALPs) include the cavity searches for galactic dark matter axions [15- 
19], searches for solar axions based on the "helioscope" method [15,20-24], the Bragg 
scattering technique [25-28] or the resonant method involving nuclear couplings [29-31], 
the polarization of light propagating through a transverse magnetic field [32-35], and 
the photon regeneration method [34, 36-45] . 

Almost all of these past, present or future experiments rely on the axion coupling 
to two photons that is described by the Lagrangian term 

^a'y = -\ g.^F^.F^^'^a = (7a7 E ■ B a , (1) 

where a is the axion field, F the electromagnetic field-strength tensor, F its dual, E the 
electric, and B the magnetic field. The axion-photon coupling strength is quantified by 
a 1 /E 2 4: + z + w\ 

^a7 = 7T- X U? - o 1 , ^ , ... ' (2) 



271 \N 3 1 + z + w. 
where z = m^/m^ and w = m^/rris <^ z are quark-mass ratios. The canonical values 
frequently used in the axion hterature are z = 0.56 and w = 0.028 [46,47], although, 
for example, z could lie in the range 0.3-0.6 [6]. However, by far the largest uncertainty 
in the relation between g^^^ and /a comes from the model-dependent ratio E/N of the 
electromagnetic and color anomalies of the axial current associated with the axion field. 
Frequently cited generic examples are the Kim-Shifman-Vainshtein-Zakharov (KSVZ) 
model [48,49] where E/N = and the Dine-Fischler-Srednicki-Zhitnitskii (DFSZ) 
model [50,51] where E/N = 8/3, but a much broader range of possibilities exists [52]. 
The KSVZ model often serves as a prototype for a hadronic model where axions do not 
couple directly to ordinary quarks and charged leptons. 
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The second parameter relevant for axion searches is the axion mass, that is related 
to the PQ symmetry breaking scale /a by 

z'/' Um^ . ,, 10'GeV 
"^a = — T- = 6 eV , 3 

1 + 2 /a /a 

where z = 0.56 was assumed. Here, m^r = 135 MeV is the pion mass and /^r ~ 92 MeV 
is its decay constant. In terms of the axion mass, the axion-photon coupling is 

_ ^ /-E _ 2 4 + z \ l + z ma , . 

This linear relationship between gg^^y and m^, defines the "axion line" for a given axion 
model in the ga^-m^ parameter spacej]. Barring fine-tuned cancellations, QCD axions 
are expected to lie near the narrow band of parameters in the ga.-y-m.^ plane that is 
defined by the corresponding 7r° properties. 

Currently, laboratory searches for axions from the Sun [21-23,26-28] are being 
extended by the CERN Axion Solar Telescope (CAST). It takes advantage of the 
expected solar axion flux that would be produced by the Primakoff effect, i.e., the 
conversion of thermal photons in the Sun into axions in the presence of the electric fields 
of the charged particles in the solar interior. CAST itself then uses the reverse process of 
coherent axion-photon conversion in a "magnetic telescope," i.e., in the bores of a 9.26 m 
long LHC dipole magnet oriented toward the Sun. CAST has recently reported the best 
experimental limit obtained so far of g^,^ < 1.16 x 10"^^° GeV~^ for ma < 0.02 eV [24]. 
This limit is comparable to the most severe constraint from the population of horizontal- 
branch (HB) stars in globular clusters that implies g^,^ < 10^^° GeV"^ [7, 53]. 

Here we report new results from the CAST experiment obtained in 2004 with 
improved conditions in all detectors. In addition, the solar axion flux is calculated from 
a modern solar model for the purpose of the CAST data interpretation. For the first 
time we obtained better sensitivity than that arising from energy-loss arguments on 
globular cluster stars. 

The result reported here provides our final limit for the CAST vacuum setup (so- 
called CAST phase I) where the sensitivity is essentially limited to ma ^0.02 eV as 
explained below. This implies that CAST has not yet reached the "axion line" in 
the ga-y-m.^ space, but rather provides the most restrictive limits on the two-photon 
vertex of axion-like particles that are somewhat lighter for a given interaction strength 
that expected for QCD axions. Meanwhile, CAST has been refurbished as a tuning 
experiment, planning to explore the mass region up to about 1 eV (CAST phase II). 
CAST will then probe QCD axion models and will also be sensitive to the existence 
of large extra dimensions [54], introduced in the world-brane scenarios to solve the 
hierarchy problem in particle physics [55]. 

We begin in section 2 with a new calculation of the solar axion fiux. In section 3 we 
briefiy describe the CAST experiment. In section 4 we report the 2004 measurements 

ft Depending on the value of E/N, the axion-photon couphng can be positive or negative. In the 
foUowing we always mean Ig^-yl when we write g^-y^ i.e., we always take gi^^ to be a positive number. 
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and their analysis. In section 5 we present the combined result of the CAST phase I 
and conclude in section 6. 

2. Solar axion flux 

2.1. General expression 

Particles that interact with photons by the Lagrangian of equation (1) are produced in 
a hot thermal plasma by virtue of the fluctuating electric and magnetic fields, i.e., the 
fluctuations of E ■ B in the plasma act as a source for the production of axions or similar 
particles. Scalar particles couple to (E^ + B^) instead of E • B for pseudoscalars. Hence 
for scalars, the fluctuations of E^ and B^ would play the analogous role. In principle, the 
axion emission rate for a given temperature, density and chemical composition can be 
calculated directly from the spectrum of E ■ B fluctuations as was done for a degenerate 
plasma in [56]. For the purpose of calculating the solar axion flux, this formalism 
gives the same result as the simpler approach where one calculates different "processes" 
contributing to the emission rate [57]. The dominant contribution derives from those 
fluctuations where E is provided by the charged particles of the medium whereas B 
comes from propagating thermal photons. This is the usual Primakoff process where a 
photon converts into an axion in the electric field of a charged particle. Of course, in 
a plasma the charged particles are correlated in such a way that charges are screened 
beyond a certain distance, an effect that is crucial to take into account. Other emission 
processes that are subdominant include the "electro-Primakoff effect" where the B field 
is provided by a moving electric charge, i.e., the axion is produced from two virtual 
photons, each attached to a charged particle [57]. In the Sun, all particles, including the 
electrons, are nonrelativistic so that the B fields associated with moving electric charges 
are small, explaining that the electro-Primakoff effect is much less important than the 
ordinary Primakoff effect involving real photons. 

Axions or similar particles typically will also interact with other particles in the 
plasma, notably the electrons and nuclei. However, in the absence of a signal it is 
conservative to use only those processes at the source that are implied by the very 
interaction structure that is used for the detection, i.e., the axion-two photon vertex. 
Moreover, for QCD axion models, the solar fluxes derived from other processes are 
expected to be much smaller than the Primakoff-induced flux if we take account of 
stellar energy-loss limits on these interactions. Therefore, our results are especially 
relevant for hadronic axions that do not couple to electrons at tree level so that the 
Primakoff process is expected to be the dominant process and a non-observation of 
axions by CAST provides new and nontrivial limits. 

The transition rate for a photon of energy E into an axion of the same energy by 
the Primakoff effect in a stellar plasma is [58] 
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where T is the temperature (natural units with h = c = = 1 are used). Recoil 
effects are neglected so that the photon and axion energies are taken to be equal. The 
screening scale in the Debye-Hiickel approximation is 



2 Aira ( 



nuclei 

where is the electron number density and the number density of the j-th ion of 
charge Zy 

The total axion number flux at the Earth (average solar distance Dq from the 
Earth) is 

where /b = (e^/"'" — 1)~^ is the Bose-Einstein distribution of the thermal photon bath 
in the solar plasma and r = R/Rq is a dimensionless solar radial variable, normalized 
to the solar radius Rq. We take into account the plasma frequency for the photons in 
the system 

"^■ = — ■ 

It enters the dispersion relation as E'^ = k"^ + cUpj so that dk/dE = E/k. Of course, the 
value of the plasma frequency depends on the radial position in the Sun. 

The approximations used in this calculation give us the axion flux for the range 
of X-ray energies relevant for CAST and also give us the total number flux and total 
energy loss rate, with an estimated precision of a few percent. However, at energies 
near and below a typical solar plasma frequency, i.e., for energies near or below 0.3 keV, 
this calculation is not appropriate because the charged particles were treated as static 
sources of electric fields, neglecting both recoil effects and collective motions. The 
charged particles of the plasma themselves undergo plasma oscillations with frequencies 
near ujpi. Therefore, even though the thermal photon energies are cut off below ujp\, 
the emitted axion energies can reach down to zero [58], although with a steeply falling 
spectrum. For applications where this low-energy flux might be of interest, perhaps 
the small solar axion flux at optical energies, a new calculation is needed that properly 
includes the dynamical aspects of the collective plasma motions, not only the static 
screening effects that are included here. 

The X-ray telescope used in CAST is an imaging device so that we do not necessarily 
want to treat the Sun as a point source. We determine the differential solar axion flux 
as an apparent surface luminosity (pa,{E, r) of the solar disk. Put another way, we 
determine the flux per unit surface area of the two-dimensional solar disk as a function 
of the dimensionless radial coordinate < r < 1. The total axion number flux at the 
Earth is 

/•I rOG 

$a = 27r / drr / dE^^{E,r), (9) 

Jo Ji^pi 
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i.e., (fa is in units of cm~^ keV~^ per unit surface area. The unit surface area is 
dimensionless because the radial coordinate r on the solar disk is dimensionless. Of 
course, instead of using a dimensionless variable r that is normalized to the solar radius, 
we could use one that is normahzed to the apparent angular radius of the Sun in the 
sky. In this case, (pa_ is to be multiplied with {Rq/DqY, giving us the surface flux per 
steradians of solid angle and therefore the dr-integration in equation (9) would extend 
to the angular radius Rq/Dq of the Sun. 

To determine the surface luminosity we re-write the integral over the solar 
volume as 

■ / dr — s> 27r / drr dz . (10) 

Jo Jo 7'"- °f ^ ^ 

sight 

We assume the Sun to be projected on a 2-dimensional disk, i.e., we neglect parallax 
effects, where the ratio Rq/Dq ^ 0.5% sets the accuracy level of this approximation. 
We may then write the volume integral as 

dz = 2 dp , f , , (11) 

line of I' _ ^-2 V J 

sight V r 

where p is the dimensionless radial position in the Sun where the physical variables 
(temperature, density, etc.) have to be taken from a solar model. Altogether we find 

where /b and Fa^-y are to be taken at the position p. 



2.2. Integration over solar models 

In order to perform these integrations we use a recent solar model for which detailed 
data are publicly available, the 2004 model of Bahcall and Pinsonneault [59], which is 
tabulated on a fine grid (1071 points). The axion flux parameters are found to be 
$a =^io3-75 X 10^^ cm~2 g-1^ 

La = glo 1.85 X IO-^Lq , 
{E) = 4.20 keV , 

(E^) = 22.7 keV^ , (13) 

where g^ = g.^^/lO"^'^ GeV'^ The maximum of the distribution is at 3.00 keV. For 
self-consistency we have used the values of the solar radius [Rq = 6.9598 x 10^° cm) 
and of the solar luminosity {Lq = 3.8418 x 10'^^ erg s~^) from [59]. 

In figure 1 we compare the differential axion flux obtained from this model with 
that from the 1982 model [60] that was used in an earlier calculation of the solar axion 
flux [20] . Even though a large number of details has changed, standard solar models and 
their neutrino predictions have remained remarkably robust. The same is true for the 
axion flux prediction that depends only mildly on the exact solar model. Even though 
details of our understanding of the Sun may change in the future, the solar neutrino 
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Figure 1. Comparison of tiie solar axion flux calculated from a modern solar 
model [59] ( ) and an older solar model published in 1982 [60] ( ). An axion- 
photon coupling of 1 X 10^^" GcV^^ is assumed. 



measurement of all active flavors [61] has provided a direct verification of the Sun's inner 
properties so that it is hard to imagine that adopting a standard solar model for the 
axion flux calculation could be a source of gross error. 

An analytic approximation to the solar axion flux spectrum may sometimes be 
useful. Instead of the analytic form proposed by van Bibber et al. [20], we find that an 
excellent fit is provided by a three-parameter function of the kind 

^*i = cf— "jV^^^^)^/^". (14) 



dE KEq. 

Here, C is a normalization constant while the fit parameter Eq coincides with the average 
energy, {E) = Eq. If we match the number flux, average energy, and width of the 
numerical spectrum with such a fit, we find (energies in keV) 

^ = 6.02 X 10^° cm-2 s-i keV~' gf, E^-^s^e-^/i'^os . (15) 
dE 

The fit is accurate at better than about the 1% level in the interval 1 to 11 keV, while 
the relative accuracy is bad at very low or very high energies. 

In figure 2 we show a contour plot of the solar surface luminosity in axions ipa{E, r) 
as a function of axion energy E (keV) and a dimensionless radial coordinate r on the 
solar disk. We also show the solar axion flux as a function of the energy for several 
values of r for an axion-photon coupling g;^^ = 1 x 10~^° GeV~^. Most of the axion flux 
emerges from the inner 20% of the solar disk. 

2.3. Do axions escape from the Sun? 

CAST can detect solar axions only if these particles actually escape from the Sun after 
production. To estimate the axion mean free path (mfp) in the Sun we note that the 
photon-axion conversion rate of equation (5) is identical (in natural units) with the 
inverse mfp of a photon with energy E based on the Primakoff process. Since we work 
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Figure 2. Left: Solar axion surface luminosity depending on energy and the radius r 
on the solar disk. The flux is given in units of axionscm^^ s^^ keV~^ per unit surface 
area on the solar disk. Also shown is the radial distribution of the axion energy loss 
rate of the Sun (dLa/di?) as well as the energy distribution of the solar axion flux 
(d<I>a/d£'). Right: Differential solar axion spectrum, derived by integrating the model 
shown on the left up to different values of r in units of the solar radius Rq. The peak 
of the spectrum moves towards lower energies if integration radius moves towards the 
outer rim of the solar disk. 

in the recoil-free approximation where the energy of the photon and axion are identical, 
equation (5) also gives us the inverse mfp for the reverse process of an axion with energy 
E to be converted to a photon, i.e., of axion absorption. As an example we consider an 
axion with energy 4 keV, near the average of the expected spectrum, and note that the 
temperature at the solar center is T ^ 1.3 keV whereas the screening scale is Kg ^ 9 keV. 
The axion mfp is then found to be Aa ~ Qio G x 10^^ cm ^ Qio^ x lO^^i?0, or about 
10~^ of the radius of the visible universe. 

Therefore, in the absence of other interactions, the axion-photon coupling would 
have to be more than 10^ times larger than the CAST limit for axions to be re-absorbed 
within the Sun. However, even in this extreme case axions are not harmless for the solar 
structure because they would then carry the bulk of the energy flux within the Sun that 
otherwise is carried by photons. Low-mass particles that are trapped in the Sun should 
interact so strongly that their mfp is smaller than that of photons [62]. Note that the 
photon mfp for the conditions near the solar center is less than 1 mm. Only particles 
with a mfp not much larger than this will leave the solar structure unaffected. They 
will cause a gross acceleration of the rate of energy transfer in the Sun and other stars if 
their mfp is much larger than this. These requirements are so extreme that for anything 
like axions the possibility of re-absorption is not a serious possibility. 

2.4- Is the CAST limit consistent with standard solar models? 

Our final limit, equation (21), on the axion-photon coupling implies that the solar 
axion flux is bounded by La < 1.3 x 10~^ Lq. Self-consistent solar models including 
axion losses were constructed in [63] . In particular, it was found that helioseismological 



An improved limit on the axion-photon coupling from the CAST experiment 



10 



measurements of the sound-speed profile do not exclude g^,^ much below 1 x 10~^ GeV~^ 
where La would be around 0.2 Lq. It was also found that for g^,^/ = 4.5 x 10^^*^ GeV~^, 
where La = 0.04 Lq, the solar neutrino flux would exceed the standard prediction 
by about 20%. While the all-flavor *B flux has been measured by SNO with a precision 
of about 8.8% [61], the uncertainties of the flux prediction are of order 20% [64]. We 
conclude that the solar axion flux corresponding to the CAST limit is far below the range 
where it would affect helioseismology or the measured solar neutrino flux in significant 
ways. For the purpose of CAST it is therefore justified to treat the solar axion losses as 
a negligible perturbation of the standard solar models. 



3. CAST experiment 

A detailed description of a proposed method of detecting solar axions by the CERN 
Axion Solar Telescope has been given elsewhere [24, 65]. Here we only recall that CAST 
uses a decommissioned Large Hadron Collider (LHC) prototype magnet with a field 
of -B = 9.0 T in the interior of two parallel pipes of length L = 9.26 m and a cross- 
sectional area A = 2 x 14.5 cm^. The magnet is mounted on a moving platform with 
low-background X-ray detectors on either end, allowing it to track the Sun about 3 
hours per day. Solar axions of a broad energy spectrum which peaks at L^ ~ 4 keV, 
essentially reflecting the thermal conditions in the solar interior, may be converted into 
real photons inside the transverse magnetic field. The probability of this conversion in 
vacuum [20] is 

p^^.=('-^X--H'4), (16) 



1 J V 2 

where q = m\/2E is the photon-axion momentum difference. Therefore, the axion 
signal should appear as an excess of photons above background in three different 
X-ray detectors, a conventional time projection chamber (TPC), a gaseous chamber 
micromegas (MM), and an X-ray telescope with a charge coupled device (CCD) as focal 
plane detector, when the Sun and magnet are aligned. A differential flux of photons of 

B \2 
.9.26 mj V9XrT / 

is expected at the end of the magnet in case of the coherent conversion which 
occurs for gL^l, providing for the CAST sensitivity to the axion mass of up to 
0.02 eV. Here d^a/dL' is the axion spectrum at the Earth, expressed by the analytic 
approximation (15). For qL^l the sensitivity is reduced owing to the photon-axion 
momentum mismatch. To search for more massive axions, coherence can be restored by 
filling the magnetic conversion region with a low-Z gas whose pressure is adjusted in 
such a way that plasma frequency, which acts as an effective mass for the photon, equals 
the axion mass [20]. As one can see from equation (17), due to the g^^ rate suppression, 
a lot of effort is required to improve the existing limits on axions significantly. 



0.088 cm-2 day-i keV"^ (7^0 ^'''''e-^/^-^o^ 
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4. Measurement and analysis 

The new data taking has been performed after dismounting and remounting the entire 
setup for the X-ray detection, with some changes in the detectors, in the materials placed 
nearby, and in the operating conditions with respect to the data taking considered in 
our previous paper [24]. 

4.I. X-ray telescope and CCD detector 

The X-ray mirror telescope with a CCD as the focal plane detector is placed behind 
one of the magnet bores at the west end of the magnet and is looking for X-rays 
from "sunrise" axions. The sensitivity of the X-ray telescope is significantly improved 
compared to the 2003 data taking period by a continuous monitoring of its pointing 
stability, which allows us to exploit the full potential of the system. As a consequence, 
the area on the CCD where the axion signal is expected could be reduced by a factor 
of 5.8 compared to 2003. In addition, after adding new internal and external shielding 
components, the background level has been reduced by a factor of 1.5 in comparison 
with the 2003 setup. The CCD detector has shown a stable performance over the entire 
2004 running period and has acquired 197 h of tracking data (when the magnet was 
pointing to the Sun) and 1890 h of background data (taken from the same area during 
the non-tracking periods). For a detailed description of the X-ray telescope design, its 
performance, and background systematics during the 2004 data taking period we refer 
to [66]. Although no significant excess signal over background could be detected with 
the X-ray telescope during the data taking period of 2004, a new upper limit on the 
axion-photon coupling could be derived. In the following sections we summarize the 
analysis techniques and results we obtained for the X-ray telescope. 

4.I.I. Expected axion signal. As described in section 2.2 the solar axion surface 
luminosity ipa,{E, r) is a function of axion energy E and the distance from the center of 
the solar disk r, see figure 2 and equation (12). The differential axion spectrum emitted 
by a circular region on the solar disk can be calculated by integrating (pa{E, r) from the 
center of the solar disk up to a radius r, resulting in spectra as shown in the right part 
of figure 2. It is apparent that not only the integral flux changes, but the peak energy 
of the spectrum slightly shifts towards lower energies for increasing r. For r = 1 the 
spectrum peaks at lowest energies. Therefore, a photon signal resulting from axion-to- 
photon conversion in the magnet bore subsequently observed by the CCD would show 
a similar energy and radial intensity dependence, after being modified by the spatial 
response function of the X-ray optics (point spread function - PSF) and the energy 
response of the CCD detector. 

To characterize the spatial intensity distribution of a potential axion signal on the 
CCD detector, we multiplied the axion solar surface luminosity by Pa^7 of equation (16) 
and folded the resulting distribution with the point spread function of the X-ray 
telescope. The corresponding intensity distribution (axion spot image) we would expect 
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Figure 3. Expected "axion" image of the Sun as it would be observed by the CCD 
detector assuming the axion surface luminosity shown in figure 2 and assuming zero 
detector background. To determine the axion spot on the CCD, the PSF of the mirror 
system and the total effective area of the X-ray telescope were taken into account. The 
count rate integrated over the region of the spot is normalized to unity. 



on the CCD chip is shown in figure 3. The axion spot image is shghtly distorted from 
a mere circular shape due to the PSF of the X-ray mirror optics [66]. In addition, the 
image is enlarged compared to an axion solar spot image multiplied with the scale factor 
of the mirror optics assuming a perfect linear imaging characteristics. The latter effect is 
a direct consequence of the finite spatial resolution of the optics, the spatial resolution 
of the detector (pixel size), and the limited pointing accuracy of the CAST tracking 
system ( <0.01°) which spread the image of a point source. For our simulations we did 
not take into account either that the PSF depends on photon energy or the fact that 
it is in general a function of the position on the detector (off-axis angle). We verified 
that both effects are of minor importance for CAST, because the image is close to the 
optical axis, i.e., we have to deal with small off-axis angles, and we expect a signal in a 
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Figure 4. Left: Signal-to- noise ratio for the X-ray telescope depending on the radius 
of the signal-spot. The axion-photon couplings vary from g^,^ = 2 x 10~^° GcV~^ 
for the lowest curve, through g^.^ = 5,9 to 16 x 10""'^'' GeV~^ for the curve with the 
largest peak value. Right: The axion flux expected in the signal-spot area on the CCD 
relative to the axion flux for an spot radius of the size of the solar disk (encircled flux). 
Two cases are shown: the encircled flux for a perfect linear optics (red line) and the 
encircled flux for a realistic X-ray optics taking into account the point spread function 
of the CAST X-ray mirror system (black line). 



relative narrow energy band (1-7 keV). 

For a fixed ma, the measured count-rate spectrum from axion-to-photon conversion 
is calculated by 

= AE' r AiE) R(E', E) 27r P ^^^^^^9.^^'') r dr dE , (18) 

Je, Jo Jo aE 

where Si is the count rate detected in detector energy channel i in units of counts s~^, 
A{E) is the effective area of the X-ray telescope in units of cm^, R{E, E') is the 
detector response function, and d$-y(£', g^-y, r) / dE is the differential spectrum of photons 
resulting from axion-to-photon conversion in the magnet for a given m^, g^^y, and the 
radius r in units of the solar radius Rq. The radius r directly corresponds to the radius 
of a circular signal area r^ we used to extract the potential signal, where r^ is given in 
detector pixel coordinates. This signal area is centered at the location where we expect 
the maximum signal intensity on the CCD chip. The effective area of the X-ray telescope 
includes the quantum efficiency of the detector, the mirror reflectivity, the influence of 
the finite size of the axion emission region on the Sun, and geometric effects due to, e.g., 
the finite diameter of the magnet bore (e.g. vignetting, see [66]). 

4- 1-2. Size of the signal extraction region. For a fixed location of the signal-spot on 
the CCD we expect that an optimum spot radius r^ exists, such that the signal-to-noise 
ratio of the CCD detector has a maximum. In order to find the best radius, we used 
the signal-to-noise ratio depending on Ts for different gi^j according to 

5" , . s(r„. a^^,) . 
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where s{rs,ga^) is the number of counts from axion-to-photon conversion expected in 
the signal region for a specific g^'y, calculated from equation (18) by summing over all 
energy bins in the 1-7 keV energy range. Further, b{rs) is the number of background 
counts we expect in a spot with the same area and location. To determine fe(rs) we used 
the background data which provides enough statistics [b{rs) ^ 200 counts] to avoid the 
optimization to be biased by low counting statistics. The background data was observed 
under the same axion sensitive conditions as during tracking while we were not pointing 
to the Sun (see right panel of figure 5). The left panel of figure 4 shows the resulting 
signal-to-noise ratio depending on r^ for different values of g^-y- We find a common 
maximum for different g^^ at a spot radius of rg = 11.5 pixels which corresponds to 
82.6% encircled solar axion flux (see right panel of figure 4). 

4- 1.3. Spectral fitting results. In order to minimize the influence of the Cu-K 
fluorescence line at ~ 8keV apparent in the background spectrum of the CCD on the 
sensitivity of the X-ray telescope (originating in materials close to the CCD chip [66]), 
we restricted our analysis to the energy range from 1 to 7keV. In total we observed 
26 counts in this energy range inside the signal-spot during axion sensitive conditions 
(exposure time 197 h). The background is deflned by the data taken from the same 
signal-spot area during the non-tracking periods, acquired under the same operating 
conditions except that the magnet was not pointing to the Sun. The spatial distribution 
of the observed events during tracking and non-tracking time is shown in flgure 5. 

The resulting low counting statistics required the use of a likelihood function in the 
minimization procedure to determine an upper limit on g^^^, rather than a x^-analysis. 
The likelihood function we used is based on a Poissonian p.d.f. and is given by 



where = 20 is the number of spectral bins, Ui the number of observed counts in bin i, 
and /ij the value of the flt function in bin i. As a flt function /Xj = Si + bi was used, where 
bi is the measured background and Si the expected number of counts in the energy bin 
i from axion-to-photon conversion. The best estimate for g^^ is obtained by minimizing 
5* = — 21n£(/i). In the large-sample limit the statistic S is x^-distributed [6]. In the 
Poissonian regime, but with relatively large, this is a reasonable assumption. The 
validity of the interpretation of S in our particular case, as well as the negligible 
influence of the statistical uncertainty of the background on the flnal result have been 
verifled with a toy Monte Carlo model. We compared the result derived with the 
maximum-likelihood estimator deflned in equation (20) with different other maximum- 
likelihood techniques based on an unbinned maximum-likelihood estimator, a folded 
extended likelihood estimator, a binned likelihood estimator that takes into account 
the spectral shape of the observed background, and an extended likelihood estimator 
with Poissonian convolution. According to our Monte Carlo simulations all methods 
are unbiased and yield identical results. We stress that in the absence of a signal 
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Figure 5. Left: Spatial distribution of events observed under axion sensitive 
conditions by the CAST X-ray telescope during the 2004 data taking period. The 
intensity is given in counts per pixel and is integrated over the tracking period of 197 h. 
Right: Background spatial distribution as observed by the CAST X-ray telescope 
during the 2004 data taking period. The intensity is given in counts per pixel and 
integrated over the full observation period of 1890 h. 

(background-dominated regime), the upper limit on g^^ does not improve when taking 
the expected intensity distribution in the axion spot shown in figure 3 into account. 
This holds as long as the effective size of the signal region is kept constant. 

The final analysis yields an upper limit on the axion-photon coupling of 
(7a7 < 8.9 X 10~^^ GeV~^ (95% CL). The observed energy spectrum together with the 
theoretically expected axion spectrum for the best fit value of (737 and for the 95% CL 
limit on g^r^ is shown in figure 6 together with the likelihood distribution of the fit. A 
summary of our results is given in table 1. 




Figure 6. Left panel: Spectral distribution (□) of the data measured with the X-ray 
telescope during the tracking exposure of 197 h together with the expectation fii for 

the best fit gaj ( ) and for the 95% CL limit on g^-y ( ) in units of total counts 

per energy bin in the signal-spot area (9.35 mm^). Right panel: Corresponding 
distribution depending on g'^^. The vertical line ( — • — ) marks the upper limit on 
gtJ95%CL). 



4.1.4. Systematic uncertainties. We studied the influence of systematic uncertainties 
on tlie best fit value of g^^ and on the upper limit on g^.^. As stated in [66] 
statistically significant variations of the background on long (months) and short time 
scales (hours to days) are not apparent from the X-ray telescope data. Even if the 
background level would be time- dependent, it would not play a significant role, since 
we measure the potential signal (signal-spot) and background (area around the signal- 
spot) simultaneously. A crucial point in the analysis is the definition of the area where 
the background is taken from. Ideally, if the background is homogeneous spatially 
distributed, the area surrounding and close to the signal spot could be used. Since we 
observed no significant difference between the background level and spatial distribution 
during tracking and non-tracking times, we used the same signal-spot area during the 
non-tracking periods to define the background. Alternatively, we used tracking data and 
different regions on the CCD outside the signal-spot area to estimate the systematic 
uncertainties due to the choice of the background definition. The overall systematic 
uncertainty is dominated by both the background definition and the pointing accuracy, 
that affects the effective area of the telescope, the location of the signal-spot and its size 
respectively. Other effects such as uncertainties in the detector calibration and magnet 
parameters are negligible. For the best fit value of g^^ we find 

(^?a\)bcstfit = (-0.31±°:i,,,, ± 0.28.,st) X 10-^°GeV-^ 

An effect of the systematic uncertainties on the upper limit on g^^^ is estimated to be 
less than ~ 5%. 



An improved limit on the axion-photon coupling from the CAST experiment 17 
Table 1. 2004 and 2003 [24] data sets included in our result. 



Run Detector Tracking Background (ga7)bostfit ± Iffstat x^in/d.o.f (95% CL) 











[LU KjC\ ) 




kjC\ ) 


2004 


CCD 
TPC 
MM 


197 
203 
196 


1890 
2616 
3000 


-o.3i±0:i 
n 21+1-2'' 


13.2/19 
17.6/17 
24.7/14 


8.9 
12.9 
12.7 


2003 


CCD 

TPC 
MM set A 
MM set B 
MM set C 


121.3 
62.7 
43.8 
11.5 
21.8 


1233.5 
719.9 
431.4 

121 

251 


0.4 ± 1.0 
-1.1 ±3.3 
-1.4 ±4.5 

2.5 ±8.8 
-9.4 ±6.5 


28.5/19 

18.1/17 

12.4/13 

6.1/13 

10.7/13 


12.3 
15.5 
16.7 
20.9 
16.7 



4.2. TPC and MM detector 

The TPC [67], covering both bores of the east end of the magnet and facing "sunset" 
axions, was operated continuously during the entire 2004 data taking period and 
acquired 203 h of tracking data and 2616 h of background data. The main difference 
between the 2003 and 2004 experimental setup were the operation of a new differential 
pumping system and an improved passive shield for the detector. The purpose of the 
differential pumping system is to decrease the effect of gas leaks towards the magnet 
and to minimize the possibility of damage to the TPC windows due to sudden pressure 
changes or break down. As a result of the smooth operation of the TPC detector, the 
exposure time during 2004 was five times higher compared to 2003. The TPC is housed 
in a low radioactive, 5 mm thick copper box inside a low radioactive multicomponent 
shield (22 cm of polyethylene, 1 mm of Cd and 2.5 cm of Pb). The entire shield is 
permanently flushed with nitrogen, creating an over pressure which decreases the radon 
contamination close to the detector. The installation of the shield has reduced the TPC 
background level in the region between 1 and 10 keV by a factor of 4.3 compared to the 
background level of 2003. The typical differential counting rate in this energy interval is 
4.1 X 10~^ counts cm~^ s~^ keV~^. One of the most important advantages of the improved 
background suppression is that the observed background is almost independent of the 
magnet pointing direction, which was not the case during the 2003 data taking period. 

The micromegas detector is placed behind the other bore at the west end of the 
magnet. The new MM [68] was specifically designed to eliminate the cross talk effects 
present at the previous detector and to improve the quality of the strips. The smooth 
operation of the improved MM detector during the 2004 data taking period, combined 
with the development of more effective off-line analysis techniques, allowed for further 
reduction of the background level by a factor of 2.5 compared to the 2003 running 
period. Due to the stability of the detector and the experiment in general, it was 
possible to accumulate 196 h of tracking data and 3000 h of background data. The 
resulting background counting rate is 5 x 10~^ counts cm~^ s"^ keV~^ in the 1-8.5 keV 
energy range. 



An improved limit on the axion-photon coupling from the CAST experiment 18 

The information of the energy scale for each detector and each run configuration is 
assured by periodical calibrations with the ^^Fe X-ray source. The hardware efficiency 
with respect to the photon energy has been simulated using the GEANT4 toolkit and 
verified by theoretical calculations and experimental data from the detectors that were 
characterized in the PANTER test facility in Munich [69]. 

4-2.1. TPC and MM data analysis and results. The results presented here were 
obtained after the analysis of the TPC and MM 2004 data sets listed in table 1. 
The analysis was performed by standard minimization. The accumulated effective 
background spectrum, properly normalized, was subtracted from the corresponding 
tracking spectrum and the resulting data was fitted with the theoretical axion signalff 
which scales with g^^ for various ma. The obtained results are given in table 1, showing 
no significant excess of axion-to-photon conversion-like events. Therefore, an upper 
limit on g.^^ at 95% confidence level was calculated for each of the data sets following 
the Bayesian scheme (shown in the last column of table 1). The corresponding fits are 
plotted in figure 7. 



xlO^ x10"^ 




Figure 7. Panels (left) and (right) show the experimental subtracted spectrum (• ) 

together with the expectation for the best fit gi^j ( ) and for the 95% CL limit 

on (7a7 ( ) for the TPC data set and MM data set, respectively. The structure 

at 3 keV in the expected spectrum reflects the change in the efficiency curves of the 
gaseous detectors due to the Ar K-edge. The vertical dashed lines indicate the fitting 
window. The fully coherent conversion regime (rria < 0.02 eV) is assumed. 

The experimental systematic uncertainties on the present limits have been studied. 
The daily variation of the background caused by the environmental radon is the 
dominant factor in the study of systematics for the MM analysis. Most of its effect is 
eliminated by using background data taken at times of the day as close as possible as the 
tracking times. The resulting upper limits for different choices of effective backgrounds 
clearly suggest that the systematic uncertainty is less than 2%. Systematic effects in 
the fitting procedure of the TPC data are considered by varying the background level 

ft The expected spectrum of axion-induced photons was calculated by equation (17) and multiplied by 
the detection efficiency curves. 
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until the null hypothesis test (in areas of the detector where no signal is expected) 
yields a result with a probability smaller than 5%. If taken as an uncertainty, this range 
corresponds to ~ 10% variation of the upper limit. 

5. Combined result of CAST phase I data sets 

The best fit values of g'^^ obtained for each of the detector's 2004 data are reported in 
table 1 together with their la error and the corresponding Xmm values and degrees of 
freedom. Each of the data set is individually compatible with the absence of any signal. 
A combined result has been obtained by multiplying the Bayesian probabilities of the 
three fits to get a combined probability function and obtaining the value of g^^ that 
encompasses 95% of its physically allowed (i.e. positive signals) part. The result of this 
operation is g^-y < 9.0 x 10~^^ GeV~^ at 95% CL. The same operation, including in the 
multiplication also the Bayesian probability functions of the 2003 data sets [24] listed 
in table 1, gives a final limit for the CAST vacuum setup of 

< 8.8 X 10"^^ GeV"^ (95% CL) . (21) 

This value is valid for ma ^0.02 eV where the expected signal is mass independent 
because the axion-photon oscillation length far exceeds the length of the magnet. For 
higher ma, the overall signal strength diminishes rapidly and the spectral shape differs 
(see equation (16)). The described procedure was repeated for different values of ma 
to get the entire 95% CL exclusion line. The region excluded in the ga^y^^a. plane is 
shown in figure 8 together with the results of other experiments as well as astrophysical 
and cosmological constraints. The statistical limit on the sensitivity is mostly set by 
the value of {BLY in the LHC test magnet. In general, the systematic uncertainties 
on the background spectrum due to its position and time dependence and due to the 
pointing accuracy are estimated to have an effect of less than ~ 5% in the final upper 
limit obtained, that is dominated by the result of the X-ray telescope. 

6. Conclusion 

We have searched for the coherent production of photons in a decommissioned LHC test 
magnet of L = 9.26 m and i? = 9 T, assuming the existence of an incident axion flux 
that is emitted from the Sun by the Primakoff process. Improvements to this experiment 
relative to the published 2003 results come from several changes in the detectors, in the 
shielding, and in the operating conditions. 

In conclusion, we obtained the currently best experimental limit on the axion- 
photon coupling which is by a factor 1.3 more restrictive than the CAST limit previously 
reported [24]. For axion masses less than 0.02 eV it supersedes the previous astrophysical 
limit based on the helium-burning lifetime of HB stars, as shown in figure 8. Our 
sensitivity does not yet touch the band of 5'a7^'^a values expected for QCD axion models 
and thus applies to somewhat lighter axion-like particles. 
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Figure 8. Exclusion plots in the axion-photon coupling versus the axion mass plane. 
The limit achieved by the present experiment (CAST phase I) is compared with 
other constraints (Lazarus et al. [21], SOLAX [26], COSME [27], DAMA [28], Tokyo 
helioscope [22,23] and HB stars [7,53]) discussed in the introduction. The vertical 
red line (HDM) is the hot dark matter limit for hadronic axions rria < 1.05 cV [77] 
inferred from observations of the cosmological large-scale structure. The yellow band 
represents typical theoretical models with \E/N — 1.95| in the range 0.07-7 while the 
green solid line corresponds to the case when E/N = is assumed. 



Our limits can be relevant, for example, in comparison with the anomalous 
signature, far in excess of the QED expectation of vacuum birefringence, that has 
been reported by the PVLAS experiment [35]. An interpretation of this signal in 
terms of axion-like particles requires a coupling strength of about 10~^-10~^ GeV~^ 
for TTia ~ 10~^ eV. However, this coupling is so large that the solar axion production 
would lead to the Sun burning out within about a 1000 years. Therefore, if the particle 
interpretation of PVLAS is correct, something about the underlying physics must be 
very different than assumed here [70-75]. Without a specific underlying theoretical 
model it is impossible to compare PVLAS and CAST in generic terms — the discrepancy 
is simply too large. 

In order to extend our search to QCD axions, our setup has already been upgraded 
and the data collection has started in late 2005. The CAST experiment, being operated 
in a scanning mode in which the He gas pressure is varied in appropriate steps, will 
allow us to explore the range of possible axion masses up to about 1 eV. For the first 
time an experiment providing a sensitivity to g.^,^ of the order of that derived from HB 
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stars will be able to enter into the region, shown in figure 8, which is favored by axion 
models. In this mass range axions would provide a hot-dark matter component of the 
universe, similar to neutrinos [76, 77]. Another interesting feature is that in this regime 
of operation, when the gas is inserted in the magnet pipes, CAST could be also sensitive 
to the existence of (two) large extra dimensions with the compactification radius down 
to around 200 nm. 
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